Martin representation and Relative Fatou 
cN ■ Theorem for fractional Laplacian with a 

^ '■ gradient perturbation 

u • 
^ ! Piotr Graczyk 

1/^ ■ LAREMA, Universite d' Angers, 2 Bd Lavoisier, 

49045 Angers Cedex 1, France 
^ , ' Piotr.Graczyk@univ-angers.fr 

<t^ i Tomasz Jakubowski 






en 

en 

o 

(N 



% 



Institute of Mathematics, Wroclaw University of Technology, Wyb. Wyspianskiego 27, 

50-370 Wroclaw, Poland 
Tomasz.Jakubowski@pwr.wroc.pl 



CN ! Tomasz Luks 

> 

^+ . LAREMA, Universite d'Angers, 2 Bd Lavoisier, 

\o • 

^•^ • 49045 Angers Cedex 1, France 



luks@math.univ-angers.fr 



Abstract 

Let L = A°/2 + b -V with q € (1,2). We prove the Martm rep- 
resentation and the Relative Fatou Theorem for non-negative singular 
L-harmonic functions on C^'^ bounded open sets.|^ 



1 Introduction and Preliminaries 

Analysis of harmonic functions related to fractional powers A"/^ of the 
Laplace operator is an important topic, intensely developed in recent years. 
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also for perturbations of the operator A"/^, see, e.g., |2l[6l|2] 
and references therein. 

From the probabihstic point of view, stable stochastic processes with gra- 
dient perturbations on M'^, d > 2, i.e. with the infinitesimal generator 

L = A"/2 + 5.v, (1) 

where a G (0,2), constitute an important class of jump processes, intensely 
studied in recent years. Their most celebrated case are the Ornstein-Uhlenbeck 
stable processes with b{x) = Ax, A G M. They have important physical and 
financial applications and form a part of Levy-driven Ornstein-Uhlenbeck pro- 
cesses, cf. [TilliT]. 

The motivations of this paper were to: 

(i) establish the theory of the Martin representation for singular L-harmonic 
non-negative functions 

(ii) study boundary limit properties of L-harmonic functions and to obtain 
a Relative Fatou Theorem for them 

(iii) develop the theory of Hardy spaces of L-harmonic functions. 

The topics (i) and (ii) are addressed in this article and the subject (iii) in 
a forthcoming paper. 

All these topics are fundamental for the knowledge of L-harmonic functions. 
The topics (i) and (ii) are well developed for fractional Laplacians. The Martin 
representation was established in this case in [TJ [201 SSI EH [12], see also |35| 
for a general setting of Markov processes. The Relative Fatou Theorem was 
proved for a-harmonic functions on C^'^ sets in [H], on Lipschitz sets in |37] . 
and on the so-called n-iat sets in [30], see also the survey |S1 Chapter 3]. 
Furthermore, some important variants of stable processes such as relativistic, 
censored and truncated stable processes were studied from the point of view 
of topics (i) and (ii), see [161 [T71 [SH [33] and [H Section 3.4]. Nevertheless, the 
methods of these extensions do not apply to the operator L of the form ([T]). 
Let us notice that all our results are also true for Ornstein-Uhlenbeck stable 
processes. 

On the other hand, Martin representations and boundary properties of 
harmonic functions were widely studied in the case of diffusion operators, see, 
e.g., [H [13 [31 [m [221 [2S1 [SI m [1] for the resuhs on the classical Laplacian 
A = ^f^^ d'^/dx'f, and [ini [211 [221 [ini [H [2S] for its various generalizations. 

Let us mention that the methods of this article give also interesting new 
results for operators different from L. In the case of Laplacians with gradient 



perturbations, i.e. a = 2, we get new perturbation formulas for the Green 
function and the Martin and Poisson kernels (Section 13. 4p . 

The potential theory of stable stochastic processes with gradient perturba- 
tions was started in the Ornstein-Uhlenbeck case by T. Jakubowski |27t |28] . 
Next, in the general context of gradient perturbations and a > 1, with a 
function b from the Kato class /C^~^ it was developed by K. Bogdan and T. 
Jakubowski [TII| [TT] and by Z.-Q. Chen, P. Kim and R. Song jTS]. Our work 
is a natural continuation of the research presented in |11) . 

In particular, we send the reader to [11] for the definitions of the fractional 
Laplacian, C^'^ sets. Green functions and Poisson kernels, with respect to both 
operators A"' ^ and L. The definitions of a-harmonic, regular a-harmonic and 
singular a-harmonic functions can be found e.g. in the monography [HI page 
61] and are analogous for L-harmonic functions. 

Throughout this paper, like in [11], we suppose 1 < a < 2, unless stated 
otherwise. We consider an open set D of class C^'^ and a vector field b G IC2~^ 



on W^ i.e. 



limsup / \b{z)\\x - z\''-'^-'^dz = 0. 

The potential theory objects related to the operator L defined in ([T]) will be 
denoted with a tilde ~, while those related to the operator A"/^ will be denoted 
without it. In particular Gd is the Green function of D for L and Gd is the 
Green function of D for A"/^. We fix throughout this paper a point xq E D 
and define the Martin kernel of D for A"/^ by 

Moix, Q) = lim -^4^' xeD,QedD. 
The L-Martin kernel is defined by 

Md{x, Q) = lim -£^(^1^, xeD,QedD 
y^QGoixo^y) 

and we show in Section |3] its existence. 

The starting point of the research contained in this paper are the following 
mutual estimates of Green functions and Poisson kernels for L and A"/^ (see 
[TTl Theorem 1 and (72)]). 

Comparability Theorem. There exists a constant G = G{a, b, D) such that 
for all x,y & D and z G {DY, 

G~'Gn{x, y) < Go{x, y) < GGd{x, y), (2) 



C~^Pd{x, z) < Pd{x, z) < CPoix, z). (3) 

One of the main elements of the proof of 02]) is the following perturbation 
formula, that will be also very useful in our present work (see [TTl Lemma 
12]). 

Perturbation formula for Green functions. Let x,y E M'^, x ^ y. We 

have 

GD{x,y) =GD{x,y) + / GD{x,z)b{z) ■VzGD{z,y)dz. (4) 

Jd 

We start our paper by proving in Section 12.11 a generalization of the Compa- 
rability Theorem: according to Lemma [H the constant C in the estimates ([2]) 
may be chosen the same for sets D^ sufficiently close to D. The same phe- 
nomenon holds also for the Poisson kernels PD{x,y) and Pu{x,y). In Section 
I2.2l we prove a uniform integrability result, that will be needed in proving the 
main results of the paper, contained in Sections [3] and IH 

In Section [3] we develop the Martin theory of L-harmonic functions. We 
prove the existence of the L-Martin kernel which is L-harmonic (Theorems [8] 
and [T2|) . Next we obtain the Martin representation of singular L-harmonic 
non-negative functions on D, see Theorem [T3l 

The formula (j4]) allows us to prove very useful perturbation formulas for 
Martin kernels flT^ . Poisson kernels f lTTl) and singular a-harmonic functions 
(l28l) . In Section 1331 (jl]) and ( TT^ are proved in the diffusion case a = 2. Also 
a perturbation formula (|33|) for the L-Poisson kernel is derived. 

Section m is devoted to an important fine boundary property of singular L- 
harmonic functions: the Relative Fatou Theorem (Theorem [23]) . We provide 
a proof of this theorem based on the perturbation formula for singular a- 
harmonic functions 



2 Preparatory results 

In this section we prove some results, interesting independently, that will 
be useful in proving the main results of the paper, coming in the next sections. 

2.1 Uniform comparability of Green functions and Pois- 
son kernels 

In what follows, R"' denotes the Euclidean space of dimension d > 2, dy 
stands for the Lebesgue measure on M"'. Without further mention we will only 
consider Borelian sets, measures and functions in M.'^. By x ■ y we denote the 



Euclidean scalar product oi x,y G M'^. Writing f ^ g we mean that there is 
a constant C > such that C~^g < f < Cg. As usual, a Ab = min(a, b) and 
aW b = max(a, b). We let B{x, r) = {y e W^ : \x — y\ < r}. For U C W- we 
denote 

Su{x) = dist(a;, [/") , 

the distance to the complement of U. 

In what follows D is a bounded C^'^ open set. 

For r > define 

Dr = {x e D : 5d{x) > r} . 

When r is sufficiently small, then Dj. is also a C^'^ open set, see |38l Lemma 
5], and one may show that the localization radius of Dj. varies continuously 
with respect to r. 

In the sequel we will often use the estimates of the Green function (|34). 
[19], see also |26]) of a C^'^ open set 

V \y-zr J 

and of the Martin kernel (|2Uj) 

Moreover, in the stable case, the estimates ([5]) are uniform when we consider 
the sets D,. sufficiently close to D, i.e. there exist constants c, eo > depending 
only on D and a such that for all r G [0, eo] and x,y E D^. we have 



see fIE\ Theorem 21] and \68\ Lemma 5]. We will now show analogous unifor- 
mity of constants for the fractional Laplacian with a gradient perturbation. 

Lemma 1. (i) There exist constants c, eo > depending only on D and a 
such that for all r G [0, eo] and x,y E Dr we have 

c'^GoXx.y) < GDrix,y) < cGoX^^y)- 



(a) There exist constants C, eg > depending only on D and a such that 
for all r G [0, eo], x G D^ and y G -D^ we have 

Proof. In order to show (i) we follow the proof of the Theorem 1 in [TT]. We 
analyse below its crucial points. 

1. Comparison of Green functions Gs{x, y) and Gs{x, y) for "small" sets S, 
[TT| Lemma 13], based on estimates from [TT| Lemma 11]. Thanks to property 
([7|), we see that the comparison of Green functions Gs^{x, y) and Gs^ix, y) for 
small sets S holds with a common constant c, when r G [0, eg]. 

2. Harnack inequalities for L and the Boundary Harnack Principle, [TT] 
Lemmas 15, 16]. Thanks to 1., we get them uniformly with respect to r G 

[0,eo]. 

3. Now the proof of (i) for any C^'^ open set D is the same as in Section 5 

of[n]. 

The part (ii) is implied by (i), applying the Ikeda-Watanabe formula for the 
Poisson kernel Pd, see [TT| Lemma 6 and (39)]. Recall that the Levy system 
for the process Xt is given by the Levy measure of the a-stable process Xt. D 

An immediate consequence of Lemma [1] and |26[ Theorem 22] is the following 
uniform estimate of the Poisson kernels of L for Dr- 

Corollary 2. There exist positive constants C, eg depending only on D, a and 
b such that for all r G [0, eo], x G -D,. and y G D^ we have 

C-'stix) C6t\x) 

< PoAx^y) < 



2.2 Derivatives of the Poisson kernel for A"/^ 

In this section we prove useful gradient estimates for the Poisson kernel of 
A"/2 for D, < «< 2. _ 

Consider a ball B = B{^q, r) G B G D and let Pg be the Poisson kernel of 
B 



PB{x,y)=C[ 



r^ — |x — ^oP"' 



-d 



X -yl ", X G B,y G B" 



and equal to elsewhere. By [131 Lemma 3.1], 

|V.PB(x,y)| < (rf + «)-^MM) , xGB,yG (By. (9) 

r- \x-^o\ 

We will now show analogous estimate for C^'^ bounded open sets. 
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Lemma 3. Suppose < a < 2. Let D he a hounded C^'^ open set in W^ and 
let P]j{x,y) he the Poisson kernel o/ A"/^ for D. Then we have 

\V.PD{x.y)\<{a + d)^^^, xeD,ye(DY. (10) 

Proof. For x E D denote B^ = B{x, Soix)). In view of [T2| (29)] we have 

PDix,y) = PB,ix,y)+ / PB^ix,z)PDiz,y)dz. 

By (jni) and bounded convergence we have 

\V,PD{x,y)\ < \V,PBAx,y)\ + |V, / PB^{x,z)PD{z,y)dz\ 



B'^ 



<{d + a)^^^+ [ \V,PBAx,z)\PD{z,y)dz<{d + a/''^'''y^ 



^oix) Jbc^ " ' ' 5d{x) 

U 

From ([9]) and the dominated convergence theorem it follows, that if / is 
a-harmonic in D then 

Q^Hx) = j ^ g^Psix, y)f{y) dy, z = l,...,d. (11) 

The estimate (Q and fITT]) gives ( |13[ Lemma 3.2]) 

Lemma 4. Let U he an arhitrary open set in M.'^ and let a G (0, 2). For every 
non-negative function u on W^ which is a-harmonic in U, we have 

\Vu{x)\ < d^^^ , xeU. (12) 

du{x) 

Since Gu{-,y) is a-harmonic in U \ {y}, for every y E U we obtain 

\V.Gu{x,y)\<d ,^,^f;^^ | , x,yeU,x^y. (13) 

du{x) A\x-y\ 



2.3 A uniform integrability result 

One of the important results of [IT] is 



Lemma 5. Gi:){y,w)/[5{w) /\ \y — u;|] is uniformly in y integrable against 
\b{w)\dw. 

In the next lemma we will show a similar property for the family of functions 

Lemma 6. Let x E D be fixed. There exists N = N{D, x) G N such that the 
functions 

are uniformly in Q E dD and n > N integrable against \h{w)\dw. 

Proof. In view of the properties of D and of the estimates of Or, _^ (x, w) and 
Md{WiQ)i we can choose N = N{D,x) G N sufficiently large, such that for 
all n > A^ we have 

Gd^_„ (x, w)Mn{w, Q) ^ cln^_^ {w)5n{wr^' 



5n „H - 5n _„(w)l-"/2|w - Q\'i\w - X 



< Cln „ (w) \ \w — X 



Sniwr/' 



where c and c depends only on D, a and x. The first term in the parentheses is 
integrable against \b{w)\dw independently of Q, n, so we only need to consider 
the second one. For w G D2~n and Q G dD we have 

/^/^^"f ^,, < diam(D)"/^2(^+i)('^+i-"/^), 

and for w G -D2-" \ D2-N , Q G dD we get 



< 2"/2 ( \w - 01°-'^-^ + = I 

Since Id _„(w^)|w^ — Q\^~'^~^ is uniformly in Q^n integrable against ^{w)\dw^ 
we can restrict our attention to the function 

Hn{w,Q) 



V„Hi-"/2k-Qh 



Let R> N, R eN. For k,m,n eN, k> R, m> N, n> N we define 
W^JQ,R) = \wED2-r.:^<SD (w)<^, ^-<|M;-g|<— I 

We note that WJ^^^Q, i?) = for fc < m or m > n. WJ^^^Q, R) can be covered 
by ci(2'^~™')''~^ balls of radii 2~'', where Ci = Ci{D). For r > denote 

Kr=supl \b{w)\\z-w\''-'^~^dw. 

Then i^^ — ^ as r | 0. We have 

if.„(M;,(5)|&(w)|c?w 

< (2*^+1) i-"/2(2"+i)'^2-""/2c^(2'=-"^)'^-i sup f \h{w)\dw 

^ (^2^'+i)i-"/2(2'"+i)^2-"°/2ci(2'=-"*)"'-^(2'=)°-°'-^i\:2-fe 

< C2i^2-«(2'')"/^^^2™2-""/2^ 

where C2 = C2{D,b,a). Furthermore, W^{Q,R) can be covered by C3(2^Y~^ 
balls of radii 2"*^, where C3 = C3{D), and thus 

< (2^-+i)i-"/22^'^2"""/2^3(2^y"'sup / |6(w)|rfw; 

^ (^2^+i)i-"/22^^2~""/2c3(2*^)'^"^(2'=)"-'^-^i^2-fc 
< C4/^2~«(2'=)°/^"^2-"°/^ 
where C4 = C4(£',6,«). Let A^ = {u; g L'2-" : '^d,_„ (ii') < 2^^}. Then 

00 n— 1 00 

fc=i? m=Ar k=RVm 



and we obtain 

Hn{w , Q)\b{w)\dw 



^S 



oo n — 1 oo 

k=R m=N k=RS/m 



n-1 



where cg = C6(-D, &, a). For w G -D2-" \ ^r we have 

H^{w, Q) < 2-/2(2^)1-/2 (^^ + i^) < 4'^^, 

so 5^(Q) := {w : i/„(w;, Q) > 4^^} C A^ for all Q e dD and n > N. There- 
fore 

lim sup / Hn{w,Q)\h{w)\dw < lim C6-ft'2-« = 0. 

■R"^°° QedD,n>N JB'^iQ) ^-^°° 

D 



3 Martin kernel and Martin representation 

In this section we will discuss first the existence and the properties of the 
Martin kernel of L for a C^'^ bounded open set D. Next we will investigate the 
Martin representation for non- negative singular L-harnionic functions on D. 

3.1 Existence and Perturbation formula for the L-Martin 
kernel 

In order to prove the existence of the L-Martin kernel, we will need the 
following property of the Green function for A°/^. 

Lemma 7. For all x E D and Q G dD we have 

hm -—— — = \/,^Md{x,Q). 

y~^Q GD{xo,y) 
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Proof. Let z ^ D,Q & dD and choose r > such that B(z,r) C D and 
B{z,r)r\B{Q,r) = 0. Since Goi-,y) is a-harnionic in B{z,r) for y G B{Q,r)n 
D, by (ITT]) , we have 



V7 f TD ( ^GD{w,y) 

n ( \ ^^ PB(z,r){x,w)-—- -dw 

= / \/ccPBiz,r)ix,Wj 

JB{z,rY '^D[Xo,y) 



dw, X G B{z,r). 



Furthermore, by (Q and ([5]), 

IV7 p / x|G'z)(u7,y) , ^PB(z,r){x,w)GD{w,y) 
\V^PB(z,r)ix,w)\-—- < G 

GD{xo,y) r-\x-z\ doiyr^ 

We now use the estimate [TT| (25)] and by considering the cases S^iw) > {w—yl 
and S£){w) < \w — y\ we get ^/ W2 ^ C\w — y\°'^'^~'^. Hence the last term i 



Soiy) 



IS 



uniformly in y G B{Q,r/2) fl D integrable against dw, and thus 

hm / V^PB{z,r){x,w)-—— rdw 

y^QjB[z,rr GD{xo,y) 



V^PB{z,r){x,w)MD{w,Q)dw = V^Md{x,Q). 

lB(z,rY 

The last equality follows from f ITT]) and the a-harmonicity of the Martin kernel. 

D 

Thanks Lemma [7] we obtain the main result of this subsection. 

Theorem 8. Let x E D and Q G dD. Let Md{x,Q) be the Martin kernel of 
A°/2 for D. Denote 

Id{x, Q) = Moix, Q)+ f Gd{x, z)h{z) ■ \/,Md{z, Q)dz 

JD 

The function Id{x, Q) is well defined for x E D and Q G dD and Id{x, Q) > 0. 
Moreover the following limit exists and equals: 

lim ^^*^^'^) - Id{.x,Q) 



y^QGoixo^y) iD{xo,Q)' 
Thus the Martin kernel of L = A"/^ + 6 ■ V for D exists and equals 

Md{x,Q) 



1'd{xo,Q) 



Md{x,Q)+ I GDix,z)b{z)-V,MD{z,Q)dz 

D 



11 



(14) 



Proof. We divide the perturbation formula (jl]) for the Green function Gd{x^ v) 
by Gd^Xq,!/) and let y —> Q. 

The exchange of liniy_^Q and J^ is justified by Lemma 11 of jTT], see the 
formula (49) in its proof. Note that by the Boundary Harnack Principle, 
GD{xo,y) ~ Goix.y) when y G i?(Q,eo), a sufficiently small ball around Q. 
We also use the estimates ([5]), (fT3|) and ([2]). 

The exchange of liniy^Q and V^ is justified by Lemma [71 Finally 

lim ^^4^ = ^d{x, Q)+ [ Gd{x, z)b{z) ■ VMoiz, Q) = Id{x, Q). 
y-^QGD{xo,y) Jd 

The strict positivity of the function l£,{x,Q) follows from (|2]), which implies 
that there exists a > such that 

lDix,Q)>aMD{x,Q)>0. (15) 

Now we consider the quotient 

GD{x,y) _ Gnix.y) Gnixo.y) Id{x,Q) 



GDixo,y) GD{xo,y)GDixo,y) Id{xo,Q)' 
when y ^f Q. D 

Directly from the definition of M£,{x,Q) and ([2]) we obtain the following 
corollary. 

Corollary 9. There is a constant c such that for all x & D and Q G dD, 

c-^Moix, Q) < Moix, Q) < cMd{x, Q). (16) 

3.2 Properties of the L-Martin kernel 

We will now study further properties of the Martin kernel of L for D. We 
start with the following useful formulas. 

Lemma 10. Consider a C^'^ open set U G U G D. 

(i) (T*erturbation formula for the Poisson kernelj For all x E U, z G {Uy 

Pu{x, z) = Pu{x, z)+ I Gu{x, w)b{w) ■ VPuiw, z)dw. (17) 

Ju 

(a) Let Q G dD. We have the following expression for the L-Poisson integral 

of the Martin kernel o/ A°/^.- 

PuMd{x,Q):= f Pu{x,y)MD{y,Q)dy 

= Md{x,Q)+ I Gu{x,z)b{z)-VMD{z,Q)dz, x G U. (18) 
Ju 
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Proof. In the following we apply the Ikeda-Watanabe formula for the Poisson 
kernels Pu and P[/. By (jlj) and Fubini's theorem, for any x G f/ and z G U^, 



1^ y I 



^^(1,06 



u\z-y\'^+'' 



Gu{x,y)+ / Gu{x,w)b{w) ■VGu{w,y)dw 



u 



dy 



Pu{x,z)+ / Gu{x,w)b{w) ■'VPu{w,z)dw. 
Ju 



For the necessary exchanges of order of integration and derivation in the last 
formula, we apply ( TT^ . ([2]), Lemma O and bounded convergence theorem. In 
order to prove (ii), we use (i) and insert the formula flT7|) in Jj^^ Pu{x, y)MD{y, Q)dy. 
We obtain 

Pu{x,y)MD{y,Q)dy = Md{x,Q) + / Gu{x,z)b{z) ■VMn{z,Q)dz 

Ju 

In the last equality the use of Fubini theorem and the exchange of J and V 
are justified by fITU]) . (jH]), Lemma [S] and bounded convergence. D 

Lemma 11. The Martin kernel Md{', ■) is jointly continuous on D x dD. 

Proof. By Theorem [8] and the continuity of Md{-, ■), it suffices to show the 
joint continuity on Z^ x dD of the function 

f{x,Q)= [ GD{x,z)b{z)-VMD{z,Q)dz. 

JD 

Let z E D. By fITT]) and the a-harmonicity of Mo{-, Q), for r > sufficiently 
small, we have 

VMDiz,Q) = V /" PBiz,r){^,w)MD{w,Q)dy 

J B(z,rY 

= / VPB{z,r){z,w)MD{w,Q)dw. 

J B{z,tY 

From ([9]) and ([6]) it follows, that V PB{z,r){.z,w)Mj:,{w,Q) is uniformly in Q 
integrable against dw. This implies that VM/5(z, ■) is continuous on dD for 
every z E D. Let now x E D and choose r > such that -B(x, r) C -D. By ([2]), 
dSD, ([12]) and (ED, for all y G B{x, r), z E D and Q E dD, we have 



C5D{zr-^ ^ C_ 

\y - z\'^~°'\z - Q\'^ ~ \y- z\'^~'^\z - QY 



G^(y,^)|VM^(.,g)|< ,.. .,,;!,_ ^|. < ,_,.,., _^|.+,-. - (19) 
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Hence, G j:,{y ^ z)\W M j:){z ^ Q)\ is uniformly in y G B{x,r) and Q G dD inte- 
grable against \b{z)\dz, which gives the continuity of /(-, ■). D 

We will now use Lemma [TO] to show L-harmonicity of M(x, Q). 

Theorem 12. For every Q G dD the Martin kernel M(x, Q) is a singular 
L-harmonic function of x on D. 

Proof. First consider a C^'^ open set U = D.,.. We note that 



Gd^x.w) = Gjj{x,w) + I Pu{x, z)GD{z,w)dz. (20) 

By ([IS]), (HZ]), ([20]) and Fubini's theorem 

PuId{x,Q) = I Pu{x,z)lD{z,Q)dz 

= / Pu{x,z)MD{z,Q)dz 

+ / Pu{x,z) / GD{z,w)b{w) ■VMn{w,Q)dwdz 
Ju" Jd 

= Md{x, Q) + / Gu{x, w)b{w) ■ VPu{w, z)dwMD{z, Q)dz 



+ 

'D 



Pu{x, z)Gd{z, w)dz 



b{w) ■VMD{w,Q)dw 



Md{x,Q)+ / Gu{x,w)b{w) ■VMD{w,Q)dw 



u 



[Gd{x, w) - Gu{x, w)]b{w) ■ VMd{w, Q)dw 



D 



= Md{x,Q)+ / GDix,w)b{w) ■\/MDiw,Q)dw = lDix,Q). 
Jd 

Thus the function Id{x,Q) is regular L-harmonic on each set U = Dr for 
r sufficiently small. By the strong Markov property, it has the mean value 
property on each open set U G U G D. D 

3.3 L-Martin representation 

The objective of this section is to prove the following Martin representation 
theorem for non-negative singular L-harmonic functions on D. 
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Theorem 13. For every non-negative finite measure v on dD the function u 
given by 

u{x)= f MD{x,Q)du{Q), (21) 

JdD 

is singular L-harmonic on D. Conversely, if u is non-negative singular L- 
harmonic on D, then there exists a unique non-negative finite measure v on 
dD verifying l[21\) . 



Proof. The L-harmonicity of the Martin integral (12T|) and the uniqueness of 
the representation follow from Theorem [T2l Lemma [TTl (l6]), ( IT6l) and Fubini 
theorem, in the same way as in the case of the Martin representation for a- 
harmonic functions in [71 proof of Theorem 1]. We will now focus on the 
existence part. By L-harmonicity of u and by (ITT]) we have for each n 



u{x) = / PD,/„{x,y)u{y)dy 



/ u{y) PD,^„ix,y)+ GD,^„ix,w)b{w)-VPD,^„iw,y)di 



1/ 
Denote 

u*n{x)= I PD,Jx,y)u{y)dy. 



dy. 



1/n 



By dlO]), [m (72)] and Lemma [5] we have 

/ / GD,/Ax,w)\b{w)\\VPD,^„{w,y)\u{y)dwdy 

<C Gd,,„ {x, w)\b{w) I /^^y dw < oo, 

where C = C{a, b, -Di/„) > 0. Hence, by Fubini theorem 

u{x) = <(x) + / GD^/„{x,w)b{w) ■ VPD^/„{w,y)u{y)dydw. 

The function m* is a-harmonic on -Di/„, so it is different iable. In order to 
justify the exchange of J and V in the last integral we fix u; G -Di/n- Then 
by (ITU]) and [TTl (72)], for e > sufficiently small and all w' G B[w,e) and 
y G -DJ/„ we have 

\V^.Po,^S^\y)u{y)\<C ""^y^ 
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where C = C{a,b,Di/n,£) > 0. Since the last term is integrable on -DJ/,„, by 
the dominated convergence we obtain 

u{x) = <(x) + / Gd.Jx, w)b{w) ■ Vul{w)dw. (22) 

JD 

We now study the sequence u*^{x) in the same way as K. Bogdan [7] in the proof 
of the existence part of the A"/^— Martin representation, with the difference 
that in our case the function u under the integral defining u* is not a-harmonic. 
Like in [3 (2.27)] we have 



1/n 1/n -•-/" 



<(^) = / PD,jS^,y)u{y)dy= I I u{y)Ad,a .^ '^",1^+^ d^dy 



Set /i„,(dO = '^d,aGD^/,S^o,0 Id^/„ \i:^y\l+o, dyd^. Lemma [D implies that 
finiM"^) = PD,/„{xo,y)u{y)dy<C Pd,^^{xo, y)u{y)dy = cu{xo) < oo 

(recall that if u was a-harmonic, then /i„(]R'^) = u{xo))- We obtain 

f Gd,, {x,i) 



The only other property of the function u intervening in the proof of the 
existence part of the A°/^— Martin representation in [7] is 

lim / u[y)dy = 

1/n 

and it also holds in our case: the L-harmonic function u is integrable on -DJ/„ 
for every n. The sequence (/i„) of simultaneously bounded finite measures with 
support contained in D is tight. We choose a subsequence /i„^ converging to 
a finite (perhaps zero) measure /i. This choice is common for all x. Without 
loss of generality, we may suppose that (n^) is a subsequence of (2~"). The 
limit measure fi satisfies 

supp(/i) C dD. 

Exactly as in the proof of the existence part of the A"/^— Martin representation 
in [7], we deduce that for all x G -D the limit 

lim<,(a;) = u*{x) 
16 



exists and 

u*ix)= f MDix,Q)dfx{Q). (23) 

JdD 

Furthermore, in view of fITT]) . for x G -Di/n and r > sufficiently small we have 
V<(x) = v/ PB{x,r){x,y)u*^{y)dy= \/PB{x,r){x,y)ul{y)dy. 

J B{x,rY J B{x,rY 

By Lemma [U and ([9]) we have 

|VP^(.,)(a:,y)<(2/)| < C^^^^^^uiy), 

and by the dominated convergence we get Vu* (x) — ?■ Vn*(x) as A; — )■ oo. We 
also have GD^/^ix^w) /^ Gd{x,w). In order to justify the passage with the 
limit under the integral sign in (122|) with n^ instead of n we observe that the 
functions Gdi,„ {x,w)b{w) ■ Vu*^(w) are uniformly integrable on D. Clearly, 
by Lemma[T]we have c~^u*(w;) < u{w) < cu*^{w), where c does not depend on 
ra, thus u*^{w) < cu*{w). By the gradient estimates we get 

GD,Jx,w)\biw)\\Vu:{w)\<Gn,,Ax,w)\biw)\- ^ ' 



and the uniform integrability follows from fl23|) . Lemma [1] and Lemma [6l 
Therefore 

u{x) = u*{x) + / GD{x,w)b{w) ■Vu*{w)dw, (24) 

Jd 

which, using ( l23l) . becomes 

m(x) = (25) 



Moix, Q)dii{Q) + / Gd{x, w)b{w) ■ V / M^(«;, g)d/i(g)ci«;. 

aD JD JdD 

By the gradient estimates and dominated convergence we also get 
v/ MDiw,Q)dfx{Q)= [ \/MDiw,Q)dfx{Q), w E D. 

JdD JdD 

Define a measure u on dD by v{dQ) = Id{xo, Q)dfi{Q). As the function Q — )■ 
Id{xo, Q) is continuous positive, the measure v is finite positive on dD. Using 

17 



Fubini theorem in (125|) and the perturbation formula for M£, from Theorem [HI 
we obtain 



u{x)= / MD{x,Q)dv[Q). 

JdD 



n 



Remark 1. We point out that the proof of Theorem\T^ is based on the pertur- 
bation formula. In fact, the methods used in J^ in order to prove the Martin 
representation theorem for singular a-harmonic functions can not be applied 
in the present case because the Green function GD{x,y) is not L-harmonic on 
D \ {x} as a function of y. 

Corollary 14. (Perturbation formula for singular L-harmonic func- 
tions j Let v{x) > be a singular L-harmonic function on D with the Martin 
representation 

v{x)= f MD{x,Q)du{Q), xeD. (26) 

JdD 

Define a singular a-harmonic function v* on D by 

v\x)= f M^(x,Q)-f^, xeD (27) 

JdD H^O)Vj 

Then the following formula holds 

v{x) = v*{x)+ I GD{x,w)h{w)-Vv*{w)dw. (28) 

Jd 

Proof. Observe that by 0151) there exists 5 > such that 

Id{xq,Q) > 5> 

for all Q E dD. Thus the measure dn(Q) = ^ ^ q) i^ finite and the function 
V* is well defined. By the unicity of the Martin representation and the formula 
( 125|) . the function v* defined by (E7|) is the same as the function v* defined by 
a limit procedure and associated to v in the proof of the Theorem [131 Hence 
the formula (p^ holds for v and v*. It is equivalent to ([28]) . D 



Corollary 15. Let v{x) > be a singular L-harmonic function on D. The 
functions v and v* are comparable: there exists c > such that for all x E D 

c-h*{x) < v{x) < cv*{x). (29) 

Proof. We use the Martin representations (p6|) . ([271) . the Corollary [9l and the 
fact that loixo, Q) > 6 > ior all Q E dD. D 



3.4 Perturbation formulas in the diffusion case 

In the present article we exploit the perturbation formulas in the case of 
the singular operator L = A"/^ + 6-V, l<a<2. In this short chapter we 
make a parenthesis and briefly discuss the case a = 2 and d > 3, corresponding 
to the diffusion operator 

on R*^. The potential theory for such diffusion generators was studied by 
Cranston and Zhao|22], and more recently by Ifra and Riahi[25], Kim and 
Song|32] and Luks|36]. Our methods allow to enrich this theory by some new 
perturbation formulas. 

We suppose that 6 G /C^ and we assume additionally that D is connected, 
i.e. it is a domain. Recall that Cranston and Zhao|22] worked under this con- 
dition and a complementary second condition |6p G IC^^^, Kim and Song|32] 
suppressed the condition on |6p and considered signed measures in the place 
of 6. 

Proposition 16. Let L = |A + 6 ■ V with 6 G /C^. Then the following 
perturbation formula for the L- Green function Go holds if x,y G W^, x ^ y. 

Goix, y) = Gd{x, y)+ f Goix, z)b{z) ■ V,Gb(z, y)dz. (30) 

Jd 

Proof. Note that by [32l Theorem 6.2], we have the estimate 

Gn{x,y)<G\x-y\^-'', x,y eR"". (31) 

The proof of the Proposition is the same as the proof of |lll Lemma 12] in the 
case 1 < a < 2, with ( l3Ti) replacing [TT| Lemma 7]. D 

Let us mention that a perturbation formula for the L-Green function was 
proposed in [25], but under a restrictive assumption of boundedness of the 
Kato norm ||6|| of b. A simpler direct proof of the estimate (13T]) without using 
the precise estimates [32| Theorem 6.2] should be available. 

Next we obtain a perturbation formula for the Martin kernel of Laplacians 
with a gradient perturbation. 

Proposition 17. Let L = |A + 6 ■ V with b E K,^. Then the following 
perturbation formula for the L- Martin kernel M^ holds if x E D and Q G dD. 



Md{x,Q) 



Id{,xq,Q) 



Md{x,Q)+ / GD{x,z)b{z)-V,MD{z,Q)dz 
•Id 



(32) 
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where Id{xq:Q) is a continuous function on dD, equal 



Id{xo,Q) = Md{xo,Q)+ / GDixo,z)b{z)-V,MD{z,Q)dz>0. 

J D 

Proof. We follow the proof of the Theorem |8] in the case a = 2. D 

The next perturbation formula concerns the L-Poisson kernel Pd{x,Q). 

Proposition 18. Let L = |A + 6 ■ V with b E }C\. Then the following 
perturbation formula for the L-Poisson kernel P^, holds if x E D and Q G dD. 

Pd{x, Q) = Poix, Q)+ f Goix, z)b{z) ■ W,Pd{z, Q)dz. (33) 

JD 

Proof. Observe that by the formula fl30p the function Gd has the same dif- 
ferentiability properties as the function Gu- In particular the inner normal 
derivative -^^(a;, Q) exists for x E D and Q E dD. It is known (see [25| page 
173]) and possible to prove by the Green formula that 

f)G 

Pd{x,Q) = ^{x,Q). 

The formula fl33|) then follows by differentiating of the formula f l30|) in the 
direction of the inner normal unit vector n. We omit the technical details. D 

Let us finish this section by some remarks. The formula Po{x,Q) = 
-q:^{x,Q) implies, like in the Laplacian case, that the L-Martin and the L- 
Poisson kernels are related by the formula 

MUx,Q) = ^p^. (34) 

Pd{xo,Q) 

On the other hand, if we insert the formula Md{x^ Q) = p^l^ X into fl33|) . we 
obtain using (132|) 

Pd{x, Q) = Pd{xo, Q)Id{xo, Q)Md{x, Q). 

Evaluating the last equation at Xq we obtain a formula for the function Id{xq, Q) 
intervening in the perturbation formula 



and another proof of the formula (IMl) . 
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4 Relative Fatou Theorem for L-harmonic func- 
tions 

We prove in this section an important boundary property of L-harnionic 
functions: the Relative Fatou Theorem. As in the preceding sections, we 
consider a nonempty bounded C^'^ open set D. Recall the Relative Fatou 
Theorem in the a-stable case. It was proved in j37] for Lipschitz sets D. 

Theorem 19. Let g and h be two non-negative singular a-harmonic functions 
on D, with Martin representations 



g{x)= f MD{x,Q)dfx^''\Q), h{x)= f MD{x,Q)dfx^''\Q), x E D. 

JdD JdD 

Then, for ji^'^^ -almost all Q G dD, 



x^Q n[x) 



where f is the density of the absolute continuous part of /i^^^ in the decom- 
position yU*^^) = fdfi'^'^^ + Ugllg with respect to the measure n^'^\ and x ^ Q 
non-tangentially. 

Our objective in this section is to prove an analogous limit property for 
non-negative singular L-harmonic functions u and v on D. 

If we denote the integral part of the perturbation formula (128|) by 

Iv*{x) = / GDix,w)b{w) ■Vv*{w)dw 
Jd 

then we have 

U = U* -\- lu* , V = V* -\- Jj,* 

where u* and v* are singular a-harmonic non-negative functions. We write 

v(x) v*(x) 1 -I- iiiM 

V* (x) 

The limit boundary behaviors of the quotients ^^ and ^^Ijrf^ will be related 



if we control the limit behavior of the quotients "*,V and "*,} . Thus we start 

^ u*(x) v*(x) 

with discussing the properties of the quotient -j4^ for a singular a-harmonic 
non-negative function h. 
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Lemma 20. Let the Martin representation h(x) = jQ^M(x,Q)dix^^\Q) for 
some non-negative finite measure ji on dD. Then, if Q ^ supp(/i *■'*'') 

liin h{x) = 

and if Q E supp(/i^''^) and a; — )■ Q non-tang entially 

lim h{x) = +00. 

Proof. The limit in the case Q ^ supp(/i^''-') follows easily from the Martin 
representation of h and the Lebesgue theorem. In the case Q E supp(/x '■'*•') we 
use the following result of Wu |47] . 

Let / be a A-harmonic function on D, corresponding via the Martin rep- 
resentation to a finite measure /i = fx^^^ on dD. If Q G supp/i, then 

liminf /(x) > 0, 



1,1 



provided x ^ Q non-tangentially. We have, on D of class C 
fix) = f P^ix^yMdy) < c I ^^/i(rfy) 

jdD JdD F ~ y\ 

= C(5z5(x)i-t / ^£^^(dy)<C6D{xy-^ [ MD{x,y)fi{dy) 
Consequently 



dD \X y\ JdD 



n[xj > 



C6d{xY-^ 
and the second part of the Lemma follows. 

D 

Lemma 21. The quotient j4^ is bounded. More exactly, there exists c > 
such that 

C-' < 1 + frr < c- 36 

h\x) 

Proof. Observe that by Corollary [T5l and the formula ( l28l) . the quotient ^*^V 
is bounded. More exactly, there exists c > such that 

V*[Xj 

As the function Id{xo,Q) is bounded, any singular a-harmonic non-negative 
function h is of the form v* for a singular L-harmonic non-negative function 

V. a 

22 



By ([2]) , if we denote 

Jh{x) = / GD{x,w)h{w) ■ Vh{w)dw 

J D 



' D 

then 

Ih{x) ~ Jh{x) 

In particular, by Lemma [21], the quotient Jh{x)/h{x) is bounded. We prove a 
much stronger property of this quotient in the following lemma. 

Lemma 22. Let h be a non-negative singular a-harmonic function on D, with 
the Martin representation h{x) = jQ^MD{x,Q)dfi^^\Q) for a finite measure 
I^W ofi QD Then, when Q e supp(/i'^'')) andx -^ Q non-tangentially, we have 

lim :M^ = 0. 

Proof. We will show that h(x]T(w) ^^ uniformly integrable in x G -D against 
the measure \b{w)\dw. Let e > 0. Since Jh{x)/h{x) is bounded it suffices to 
show that there is 5 > such that 

GD{x,w)h{w) 

—---—^b{w)dw<e, 37 

provided \{F) < 5. Here, A denotes the Lebesgue measure on W^. First, we 
note that 

GD{x,w)h{w) 

p h{x)dD{w) 

Gd{x,w)Md{w,Q) (h)fn\\hf AU 

dfi^ >{Q)\b{w)\dw 

FJdD h[x)dD[w) 

Md{x,Q) f f Gd{x,w)Md{w,Q). 




\b{w)\dw)dfj.^''\Q). (38) 



dD K^) \Jf Md{x,Q)5d{w) 

The function q'TJ \s(w) ^^ uniformly integrable in x, y G -D against \b{w)\dw 
(see the proof of [TT| Lemma 11]). Hence, there exists 5 > such that for 
\(F) < S, 

Jp GD{x,y)dD{w) 
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and consequently 



Gd{x,w)Md{w,Q) f GD{x,w)GD{w,y) 

\b{w)\dw = / lim -— — — - — \b{w)\dw 



p Md{x,Q)6d{w) JpOBy^Q GD{x,y)5Diw 

f GD{x,w)GD{w,y) 
= hm / ——-^—\b{w)\dw<e. 

D3y^Qjp GD{x,y)dD{w) 

Now, f l371) follows from f l38l) and Martin representation of h. For Q E suppyU *■'*'', 
liraD^x-^Q h{x) = oo from the Lemma [201 Hence, by uniform integrability, 

r \Jh{x)\ ^ y f GD{x,w)h{w) 

hm , , - < c hm / ,, ., , . \b{w)\dw 



f GD{x,w)h{w) 

;/ hm —————^—\b{w)\dw = 0. 



D 



Now, we return to the Relative Fatou Theorem for L-harmonic functions. 
Let u and v be two non- negative singular L-harmonic functions on D. By 
Theorem [131 they have a Martin representation 

u{x)= MD{x,Q)dfi{Q), v{x)= MD{x,Q)du{Q), xeD 

JdD JdD 

where /i and u are two Borel finite measures concentrated on dD. 

We decompose the measure /.i into its absolutely continuous and singular 
parts with respect to the measure u 

dfi= f du + dfising 

with a non-negative function / G L^{i^) and z/(supp(/ismg)) = 0. 

Theorem 23. (Relative Fatou Theorem^ Foru-almost every point Q G dD 
we have 

lim ^ = /(g) (39) 

when X — 7- Q non-tangentially. 

Proof. We will use the Relative Fatou Theorem for the singular a-harmonic 
functions u* and v* defined according to f l271) . 

Let Q G supp(i^) \ supp(/i). Then, if x — )■ Q, v*{x) — )■ oo and u*{x) — )■ 0, 
so lim2;_j,Q ^jT^ = 0. The formulas ( I36l) and ( l35l) imply that in this case 

hm"^^) - 



v{x) 
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Let us consider the case Q G supp(z/) fl supp(/i). As 

d^^{Q) _ r d^{Q) dfisingiQ) 



Id{xo,Q) Id{xo,Q) Id{xo,Q)' 

the Relative Fatou Theorem for the singular a-harnionic functions u* and v* 
says that for z/-alniost every point Q G dD 

lim ^ = fiQ) 

x->Q V*{Xj 

when X — )■ Q non-tangentially. The formula ( 139|) then follows by the formula 
(1351) and the Lemma [22l D 
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